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Abstract — This paper considers an amplify-and-forward re- 
lay network with fading states. Amplify-and-forward scheme 
(along with its variations) is the core mechanism for enabling 
cooperative communication in wireless networks, and hence 
understanding the network stability region under amplify-and- 
forward scheme is very important. However, in a relay network 
employing amplify-and-forward, the interaction between nodes 
is described in terms of real-valued "packets" (signals) instead 
of discrete packets (bits). This restrains the relay nodes from re- 
encoding the packets at desired rates. Hence, the stability analysis 
for relay networks employing amplify-and-forward scheme is 
by no means a straightforward extension of that in packet- 
based networks. In this paper, the stability region of a four-node 
relay network is characterized, and a simple throughput optimal 
algorithm with joint scheduling and rate allocation is proposed. 

I. Introduction 

Relaying is central to wireless mesh and ad hoc networks, 
and is a potential enhancement for existing cellular networks. 
There is no "wireless network" without multi-hop communi- 
cation, and thus, understanding the role and impact of different 
relaying techniques on wireless networks is of critical impor- 
tance. There are multiple possible forwarding strategies that 
could be used by each relay node in the network, with the most 
popular ones being decode-and-forward, amplify-and-forward 
and quantize-and-forward. The forwarding strategy used at 
the relays and the resulting network capacity are inherently 
coupled, and it is not obvious which of the above forwarding 
strategies, if any, is information theoretically optimal for a 
particular network topology. For a few settings, amplify-and- 
forward has been shown to be better than decode-and-forward 
[1], and for a few others, partial decode-and-forward has been 
shown to be optimal [2]. 

Other than purely information-theoretic rate calculations, 
there are practical reasons for the use of amplify/quantize- 
and-forward strategy. Decoding a packet at the relay requires 
an entire receive chain (demodulation and decoding) to be im- 
plemented at the relay, along with an entire transmit chain (re- 
encode, remodulate and retransmit). An amplify/quantize-and- 
forward relay can bypass this chain, processing the received 
signal directly to obtain the relay output. This simplification 
greatly impacts the cost, energy usage and size of the relays, 
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and therefore, it is not surprising that the majority of the 
relays used in practice today are based on amplify-and- 
forward strategy. Also, in practice, amplify-and-forward is at 
par with decode-and-forward strategy with reduced complexity 
in processing [3]. Therefore, it is imperative to understand 
the stability region of relay networks employing amplify-and- 
forward scheme. 

For non-cooperative relaying, there is a vast body of lit- 
erature on optimal resource allocation in wireless networks 
[4]— [6] (and references therein). For cooperative relaying [7], 
[8], the stability region of relay networks using decode- 
and-forward transmission has been studied in [9], [10]. In 
the decode-and-forward paradigm, the queues still consist 
of conventional decoded data packets. In an amplify-and- 
forward relay network, each relay, however, observes an input 
sequence, whose "rate" is dependent on the state it was en- 
coded for, which is not a conventional packet. There are three 
major differences between the amplify-and-forward relaying 
and conventional node-forwarding: 

1) Each node in the network stores and ultimately forwards 
a real valued "packet". This is accomplished by quan- 
tizing the received analog signal to some finite precision 
and storing it in a buffer. This stored signal is later 
retransmitted without further processing. 

2) The "rate" or effective bits per symbol of each real- 
valued "packet" is not the same, but variable depending 
on the encoding. 

3) A high-rate real-valued "packet" can only be forwarded 
when the channel can support the corresponding rate, 
i.e., the channel state should match the "packet" rate, 
which is not required in the decode-and-forward scheme. 

We consider a four node relay network in this paper, and 
our main contributions include: 

1) We introduce a new queue-architecture to store the input 
sequences, which allows us to optimally exploit different 
states to do amplify-and-forward. We characterize the 
stability region of the four-node relay network under 
the amplify-and-forward, which in general is larger than 
that under the conventional amplify-and-forward scheme 
which does not allow signal buffering. 

2) We also propose a throughput-optimal algorithm, which 



achieves the stability region. An interesting and impor- 
tant property of our stabilizing algorithm is that it does 
not require the knowledge of the underlying distribution 
of the fading states. 

The rest of this paper is organized as follows. The next 
section introduces the system model. Section UTTI provides an 
information-theoretic achievable rate for amplify-and-forward 
relaying. Section [IV] presents the algorithm for throughput- 
optimal stable cooperative relaying over wireless networks for 
the four node example. The paper concludes with Section IVl 
The proofs are stated in the Appendix. 

II. System Model 

Our system model consists of a source, a destination and 
two relays as shown in Figure Q] The relays 1 and 2 assist the 
source s in transmitting to the destination d through amplify- 
and-forward relaying. There is no direct link between the 
source and the destination. In addition, there is link activation 
constraint such that the links from the source to the relays and 
the links from the relays to the destination cannot be active 
at the same time. This constraint arises from the following 
two system limitations: (i) a relay cannot receive and transmit 
simultaneously due to hardware limitations, and (ii) one of 
the relays (say relay 1) cannot receive from the source at the 
same time as the other relay (relay 2) is transmitting to the 
destination. In essence, we assume that the transmissions of 
the two relays heavily interfere with each other. We consider 
a discrete time model for data transmission over the links. All 
the links undergo slow fading. Further, we make the following 
assumptions: 

(1) Assumption 1: There is coordination between the relays, 
i.e., they transmit signals synchronously. 

(2) Assumption 2: The fading process is independent and 
identically distributed (i.i.d.) block fading with block 
length of T symbols. Without loss of generality, we 
assume T = 1. 

(3) Assumption 3: There is no power control over states. We 
consider an average power constraint of P per block per 
node in the network, and additive Gaussian noise of unit 
variance at each receiver in the network. 

To simplify notation, we denote the relays by n, n G {1, 2}, 
the wireless link between the source and the relay n by l sn , 
and the wireless link between the relay n and the destination 
by l n d- Let the fading state of link l sn be F sn , and the fading 
state of link l n d be F n d- We assume that F sn and F n d are 
discrete non-negative random variables which take values from 
T . In every block, let the probability of the fading state f = 

[/si, fs2, fid, hd] 6 -F* be 7Tf . 

Consider a given fading state f . If the source transmits, then 
the received signals at the relays are 



Relay 1 



Vn = VUnXs + W n , for fl = 1, 2, 



(1) 




Relay 2 

Fig. 1. A Simple Cooperative Relay Network 

transmit, then the received signal at the destination is 

2 

Vd = X! Vfnd x n + W di ( 2 ) 
n=l 

where x n denotes the symbol transmitted from the relay n, 
and yd denotes the symbol received at the destination. Here, 
w n and Wd are i.i.d. zero-mean additive Gaussian noise of unit 
variance at the relay n and the destination. 

III. Amplify-and-Forward Achievable Rate 

In the static case without any link activation constraint, 
amplify-and-forward commonly refers to the relaying scheme 
at the relays that transmit (in every time slot) scaled versions 
of the received signals in the previous time slot. The scaling 
parameters at the relays are determined as a function of the 
signal-to-noise ratios (SNRs). These scaling parameters should 
also ensure that the average power constraints at the relays 
are satisfied. For the setting explained in the previous section, 
we look at an amplify-and-forward scheme (denoted by AF 
scheme) in which the relays can transmit any of the previously 
received signal vectors or choose not to transmit. We assume 
that any received signal vector at the relays can be transmitted 
to the destination only once. 

We say that a symbol x s is transmitted by the source 
to the destination over state g = [g s i, g S 2, gid, <?2d] £ J~ 4 , 
if the source transmits during a fading state of the form 
[g s i, <7s2, *, *] and the relays transmit during a fading state 
of the form [*,*,gid,g2d\- Consider a symbol x s transmitted 
by the source to the destination over some state g. Let the 
average power used at the source be _P s g and at the relay n 
be P*. These parameters are later optimized for the state g. 
From (Q]i, (0, the received symbol at the destination is 



yd = J2 



gndPn 

g sn Pf + l 



+ w d , (3) 



where a; s has zero mean and variance Pf. From (01, it is 
straightforward to see that the maximum rate we can obtain is 

, 2 \ 



r„ = max C 



Pf {y/9s\Cl + y/gs2C 2 ) 



(4) 



where x s denotes the symbol transmitted from the source, and 
y n denotes the symbol received at the relay n. If the relays 



ci + c 2 + 1 J 
where c„ = {g n dPf)/{g S nPf + 1) and C(x) = \ log 2 (l+a:) 



Remark 1: The rate in (0]i is equal to the maximum achievable 
rate using amplify-and-forward scheme in the static case with 
fixed channel state g and full-duplex constraint [8]. 
Remark 2: The power optimization in (0]i will result in uti- 
lizing maximum power at the source and one of the relays. In 
general, the other relay will result in using lower power than 
the maximum available power [8]. 

For handling variable rate allocation, we maintain separate 
virtual queues at each relay based on the possible rates of 
the real-valued "packet". This is necessitated by the fact that 
encoding and decoding in amplify-and-forward relaying is an 
end-to-end process. Thus, variable rate coding at the source 
depending on the state must be done to ensure that: (i) the 
virtual queues at the relays are stabilized, and (ii) the relays 
possess sufficient "packets" at variable rates for throughput 
optimal operation. We denote by Qf t [t] the length of the virtual 
queue maintained for state g, where g = [g s i,g S 2, gid, gzd] G 
T A . Note that when the system is in fading state f, for any 
state g such that g s \ = f s i and g S 2 = f S 2, a "packet" 
generated with rate r g can be successfully transmitted from 
the source to the relays; and for any state g such that 
gid = fid and g 2d = hd, a "packet" generated with rate 
r g can be successfully transmitted from the relays to the 
destination. We define I s = {(f,g)|fi> s i = f s i,g S 2 = /s2>, 
and I d = {(f,g)|sid = f\d,92d = /ad}. 

The throughput region of the network is defined in the 
following lemma. 

Lemma 1: A rate r is supportable in the four-node relay 
network using an AF scheme only if there exist af > and 
bf > such that 

r =E E ( 7r f a f r g 1 {(f.g)e/ s }) . ( 5 ) 
g fE^ 4 f e.^ 4 

51 7r f a f 1 {(f,s)e^} = £ n* 6 ? 1 {(f*)e/d}> V §> < 6 > 

2 af + 6f <1, Vf, (7) 

ge^ 4 

where 7Tf is the probability that the system is in fading state 
f and denotes the indicator function of event E. When 
the system is in fading state f , the source can transmit one 
"packet" to any one virtual queue g for which l{(f,g)e/,} = 1; 
or the relays can transmit one "packet" in any one virtual 
queue g to the destination for which l{(f,g)ei" <i } = L 
Proof: Let af be the fraction of time the source transmits 
to virtual queue g when the system is in fading state f , and 
bf be the fraction of time the relays transmit the "packets" 
in virtual queue g to the destination when the system is in 
fading state f. Hence, Inequality (0 says that the overall 
fraction should be no more than one. Equality (0 is the flow 
conservation constrain for the source to all the virtual queues, 
and Equality © is the flow conservation constrain for each 
individual virtual queue. The detailed proof is given in the 
Appendix. ■ 
We can obtain rates strictly greater than the average of 
rates over all fading states by asynchronously combining states 




0.2 0.4 0.6 0.8 1 
Outage Probability (y) 



Fig. 2. Amplify-and-Forward: Asynchronous vs Synchronous 



between the source and the relays and the relays and the 
destination. We will demonstrate this using a simple example. 
Let P = 1, T = {0, 1, 10}. Consider the fading distribution 
such that fading states [0, 0, 0, 0], [0, 0, 10, 10], [1, 1, 0, 0] and 
[1,1,10,10] occur with probabilities r ) 2 , 77, 77, and j 2 , 
respectively. Here, 7 = 1 — 7 and 7 can be considered as the 
probability of outage in this example. The rate corresponding 
to state [1,1,10,10] alone is non-zero, which is C(20/ll). 
In this example, it is easy to observe that we can achieve 
(O.57 2 + 77)C(20/11) bits/transmission with asynchronous 
combining of states whereas 0.57 2 C(20/11) bits/transmission 
is the average of rates over different fading states. These 
achievable rates are plotted in Figure [2] 

IV. Throughput Optimal Algorithm and Stability 
Analysis 

We consider an i.i.d. arrival process A[t] for the data bits 
at the source s with mean rate A and bounded variance. Let 
the queue at the source be Q s with queue length Q s [t] at time 
t. At the relay n, we assume different virtual queues 
with queue length Q^[t], which is the virtual queue length 
corresponding to state g at time t. The queue Q s consists 
of bits whereas queue Qf t consists of real-valued "packets" 
(quantized to required precision) encoded for state g at rate 
r g . We will show that this queue-architecture at the relays can 
be utilized to obtain a throughput-optimal algorithm. 

In this section, we provide an algorithm, based on back- 
pressure [4] which is throughput-optimal, and that does not 
require the knowledge of the fading distribution. The algo- 
rithm has similarities with the maximum differential backlog 
algorithm for conventional networks and its generalization to 
cooperative relaying with decode-and-forward in [10]. How- 
ever, the fact that "packets" at the relays have variable rate 
introduces different weighting factors for different "packets". 
Back-pressure-based amplify-and-forward algorithm: Sup- 
pose the system is in fading state f [t] at time slot t, then the 
central coordinator computes 

81 e argmax ^Q s [t] - r s ^ QSMJ r g 1 {(f,«)e/.} 



and 



If 



g^ G argmax r g ^ Q*[i] r g l {(fjg)e/d 



- r sl J2 Qn [t}j r gl > |r g j £ Q« 3 [t]J r g . , 

then the system operates the links from the source to 
the relays and transmits r g * bits, encoded into one real- 
valued "packet", from queue Q s [t] to queues Qf 1 [t] and 
Q2 1 [t] ■ Otherwise, the system operates the links from the 
relays to the destination and transmits one real-valued 
"packet" from queues Qf 2 [t] and Q^ 2 [t] to the destination. 

■ 

Theorem 1: The back-pressure-based amplify-and-forward al- 
gorithm stochastically stabilizes the queues for any A if there 
exists e > such that A + e is within the throughput region 
given in Lemma Q] Hence, the throughput region in Lemma Q] 
is also the stability region of the four-node relay network. 
Proof: The proof is provided in the Appendix. ■ 
In the above algorithm, the weight associated with the 
virtual queue for the "packets" transmitted from the source to 
the relays in state g is (Q s [t\ - r g ^=i Q|M) 1 {(f,g)e/ 8 }' 

The sum J2n=i QfXA appears as the two relays cooperate to 
transmit the "packets". This sum has a normalization factor 
of r g as each "packet" corresponds to r g effective bits. The 
term l{(f, g )g/ 3 } arises as, given a fading state, only "packets" 
corresponding to certain states can be transmitted from the 
source to the relays. Similar comments apply to the weight 
associated with the virtual queue for the "packets" transmitted 
from the relays to the destination in state g. 

V. Conclusion 

In this paper, we characterize the maximum stable through- 
put for a two-hop cooperative relay network. A key feature 
of this paper is that we bring together physical and network 
layer constraints in characterizing this throughput. We believe 
that the analysis conducted in this work and the algorithm can 
be generalized to other cooperative relay networks employing 
an amplify-and-forward strategy. We also believe that this 
framework can be used to analyze other forwarding strate- 
gies including partial-decode-and-forward and quantize-and- 
forward. 
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Appendix 

A. Proof of Lemma [7] 

Let af be the fraction of time the source transmits "packets" 
corresponding to state g when the system is in fading state f , 
and of be the fraction of time the relays transmit the "packets" 
corresponding to state g when the system is in fading state f . 
It is clear that the maximum achievable rate r max using the 
AF scheme is 

max min \ ^ a f r s^ TTfbfr s \ , (8) 

s.t. ( a t + b f) ^ !> Vf > 

af = 0, V(f,g)£/ 8 , 
6? = 0, V(f,g)£l d , 
af,of>0, Vf,g. 

Now, we will prove that for any r < r max there exists 
af > and of > such that ©, ©, and © are 
satisfied. Let af,6f be an optimal solution to the problem 
©. Consider the following assignment. For each state g, 
let af = (r/r m ax)O s af, and of = (r / r max )r/ s bf . Choose 
9 S and r/ s such that X)fG^ 4 7r i a f r s = er* ^fb^r^ = 

(r/r max ) min |l]fG^ 4 Tf^f r g) DfgjM T&f r g|- Note that 
< 9 S < 1 and < r) g < 1. Therefore, this assignment 
satisfies (0, ©, and This completes the proof. 

B. Proof of Theorem Q] 

Consider the following optimization problem: 



max 
of,5f 



S.t. 



fe:F 4 



n* ?^. 
Y ntafr s < Y^ 

ge^ 4 

of = 0, V(f,g)£l s , 
6f = 0, V(f, S )t h, 
af,6f>0, Vf,g. 



(9) 



f ' g> 



Vg, 



a? + «)<!, Vf, 



Let the optimal objective value of problem (O be f max . We 
will show that f max = r max . Since any feasible assignment 
for (0 is a feasible assignment for (0 and has the same 
objective value, f max < r max . Since any feasible assignment 
for (0 lead to another feasible assignment for (0 with same 
objective value (as shown in the the proof of Lemma [TJ, 
r max < r max . Hence, f max = r max . The structure of this 
alternate characterization is used later in the proof. 

When the fading state at time t is f , the algorithm in Section 
HVl is based on the optimization problem: 



max Q' W - r s E Q' W )r e a$ + 

Qf f ge^ 4 V 



(10) 



n=l 



E h- g E^w r «#> 

ggJT4 \ „=1 / 

s.t. ^(af + /3f)<l, 
ge^ 4 

af = 0, V(f,g)£I s , 
/f = 0, V(f,g)£I d , 
af,/3 f s e{0,l}, Vg. 

Note that the variables are restricted to integer values {0, 1}. 

Since the queues form a Markov chain, we can use Foster- 
Lyapunov theorem (see Proposition 5.3 in [11]) in order to 
prove the stability. We assume that r g > for all g. Otherwise, 
the queues at the relays corresponding to zero rates can be 
removed without affecting the rates supportable by the system 
and the stability of the system. Consider the Lyapunov function 



V(Q[t]) = Q 2 s [t] + J2 E^M) 2 ' 

where Q[t] denotes the vector of all queue lengths. Let an 
optimal assignment to problem dTOb be af,$f and E be the 
event Q s [t] + A[t] > J2 g er 4 r g<*f ■ We nave 



Q 2 s [t + 1] = ( Q.[t] + A[t] - [ r ^ 

2 



< \Q,[t]+A[t]- J2 
ge^ 4 



r g df 



< Q 2 s [t]+A 2 [t]+ [ r s&t 
kge^ 4 



E 



g"f 



A[t] 



Similarly, we have 
(r s QS[t+l\) 2 < (r s QS[t]+r s af-rJf 



= (r s Qm) 2 + r 2 JPf-&f 



2r s QS[t]r s ($f-af 



Let af,bf be any feasible assignment to the optimization 
problem (|9). Using the law of iterated expectations (condition- 
ing on the fading state), we obtain 

E[V(Q[t + l])-V(Q[t})\Q[t}] 



< E^ 

fe^ 4 

2 



-2Q.[t] ^ r g df-A 



53 53 (2r K Q«[t]r g (#-a?))+M 

,1 = 1 ge JT4 



= 2 



a- 53 (w g fif) 



f,ge^ 4 



E E E 



M 



< 2 



Q,[t] A - 53(7r f r g a?) 



f.g 



EE* 

n=l gS^ 4 



53 (7Tfr g (a: 



&?)) 



M, 



where M is a finite value, as the variance associated with 
the arrival process is bounded and all r g are finite. The 
last inequality holds due to the following reason: Consider 
the linear program (LP) obtained by relaxing the integer 
constraints of the optimization problem dTOb . This relaxation 
is tight as LPs have at least one optimal solution which is 
a boundary point. The feasible assignment set af,bf to the 
optimization problem (O is a subset of the feasible set for the 
LP. 

We will now show that for A < r max , there is strict negative 
drift on the set of all possible queue states, except on a 
compact subset. For a given state g, if S(f g )e/ d 7r f = 0' 
then Qf t [t] = for all t, n as the algorithm will never 
choose to transmit to these queues. Let <j> = - A . Since 
A < r max , <f> < 1. Let af,bf be an optimal solution to 
the problem (0. Note that the following is another trivial 



= 0, 



feasible assignment for the problem (|9): 
for (f , g) £ Id, bf = for (f , g) ^ Id- Since the feasible set is 
convex, any convex combination is also feasible. In particular 



consider the following convex combination: 



ft 



y f — 2 W ~i 

is clear that (A 



(1-0) 



(0+1) AS 
2 "f ' 



I l{(f ;g ) 6 / d }. For this assignment, it 
J2f ge ^4 7r f r g d s ) < 0, and, for any given 

state g, Efe^Ofglof - Pf)) < if £(f, g )ei d > °- 
This completes the proof. 



